第 39 卷 第 9 期 


py 
中 


计算 机 应 用 研究 Vl 


月 十 | 
No.9 


Application Research of Computers Accepted Paper 


无 迹 西 格 玛 点 引导 的 拟 反 向 忒 菌 算法 及 其 工程 应 用 
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摘 要 : 针对 秀 菌 算法 搜索 停滞 和 算法 稳定 性 差 等 问题 ， 提 出 了 无 迹 西格玛 点 引导 的 拟 反 向 针 菌 算法 。 首 先 ， 使 用 
了 拟 反 向 学 习 和 拟 反 射 学 习 两 种 反 向 学 习 过 程 ， 根 据 原 始 和 菌 算法 勘探 和 开采 行为 的 表现 时 机 ， 生 成 同时 包含 拟 反 


向 和 拟 反射 的 


反 向 种 群 ， 扩 大 搜索 范围 ; 其 次 ， 根 据 种 群 的 多 样 性 程度 判断 是 否 使 用 反 向 种 群 重 构 原 始 种 群 进 


行 后 续 计 算 ， 避 免 固定 的 反 向 过 程 破坏 种 群 本 身 的 搜索 特点 ， 提 高 搜索 精度 ; 最 后 ， 引 入 无 迹 变 换 的 西格玛 点 ， 改 
进 黏 菌 算法 的 基本 移动 模式 , 使 无 迹 西格玛 点 引导 黏 菌 算法 的 搜索 ， 加 快 收 争 速度 。 实 验 部 分 使 用 了 CEC2017 基准 
测试 函数 ， 使 用 了 传统 统计 特征 和 MAE 排名 、Wilcoxon 秩 和 检验 验证 算法 有 效 性 ， 并 使 用 其 求解 轿车 侧面 碰撞 的 


实际 工程 优化 问题 ， 
策略 有 效 且 各 策略 间 组 合 相 


与 新 颖 的 高 水 平 群 智能 算法 、 改 进 算 法 、 不 完全 算法 进行 对 比 测试 ， 实 验 结果 表明 提出 的 改进 
改进 后 算法 的 求解 精度 和 鲁 棒 性 更 具 竞 争 力 。 
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Unscented sigma point guided quasi-opposite slime mould algorithm and 
its application in engineering problem 


Liu Yusong, Liu Sheng+ 


(School of Management Studies, Shanghai University of Engineering Science, Shanghai 201600, China) 


Abstract: Aiming at the search stagnation and poor stability of the Slime Mould algorithm, propose an unscented sigma point 


guided quasi-opposite Slime Mould Algorithm. Firstly, use quasi-opposite learning and quasi-reflected learning to exploration 


and exploitation behaviors according to the original Slime Mould Algorithm to generate a comprehensive opposite population 


that includes both quasi-opposite learning and quasi-reflected learning, expand the search space. Secondly, according to the 


diversity of the population, decide whether to use the opposite population to regenerate a new population for subsequent 


calculations, avoid the continuous opposite process destroying the search characteristics of the population itself, and improve 


the search accuracy. Finally, use unscented transformation Sigma point to improve the basic movement mode of Slime Mould 


Algorithm, make the unscented sigma point guide the search, and accelerate the convergence speed. The experimental part 


uses the CEC2017 benchmark test functions, uses traditional statistical index and MAE ranking, Wilcoxon rank-sum test to 


verify the effectiveness of the algorithm, and use it to solve the car side impact design problem, compares and tests with the 


novel swarm intelligence algorithms, improved algorithms and incomplete algorithm. The experimental results show that the 


strategies ls effective and combinations of strategies complement each other and the improved algorithm's solution accuracy 


and robustness are more competitive. 
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0 引言 
在 大 型 复杂 优 


几 年 间 愈 发 流行 


度 低 等 问题 ， 虽 然 现 阶段 对 SMA 算法 的 研究 相对 较 少 ， 
内 外 学 者 依旧 使 用 不 同 的 策略 进行 改进 获得 了 不 同 的 成 果 。 


区 
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在 大 数据 和 人 工 智能 的 背景 下 ， 群 智能 启发 式 算法 因 其 ”在 使 用 反 向 学 习 框 架 改进 黏 菌 算法 的 方向 


， 文 献 [6] 使 用 混 


ph 优秀 的 表现 和 低廉 的 算 力 要 求 在 近 十 ” 沌 反 向 学 习 框 架 扩 大 算法 搜索 空间 ; 文献 [7] 使 


] 
(slime mould algorithm, SMA) 习 框 架 帮 助 算法 跳出 局 部 最 优 : 文献 [8] 使 用 精英 反 向 和 二 次 


站 


求解 随机 优化 问题 。 插值 的 综合 策略 改进 了 笑 菌 算法 ， 使 其 整体 表现 更 优 。 在 改 


动态 自 适应 、 反 进 厅 菌 算法 基本 移动 方式 的 方向 上 ， 文 献 [3] 使 用 名 为 “ 况 食 - 
着 空间 ， 受 到 了 ”进攻 ”的 策略 , 基于 贪 梦 策 略 和 高 斯 变异 , 针对 单个 个 体位 置 


启发 而 设计 ， 通 过 生物 


菌 的 振荡 模式 ”每 轮 欠 代 中 进行 新 位 置 的 尝试 ， 使 得 求解 精度 大 大 提高 ， 文 


位 置 食物 信息 的 正 向 献 [9] 将 人 工蜂 群 的 交叉 算 子 引入 黏 菌 算 法 ,加快 算 法 收敛 速 
前 该 算法 已 成 度 ; 文献 [10] 将 正 余 弦 算 法 的 搜索 方式 添加 进 舟 菌 算法 , 提高 
下 第 选中 \、 医学 Xx 算法 的 搜索 效率 ; 文献 [11] 使 用 量子 旋转 门 替 换 黏 菌 算 法 的 


本 


实际 问 


题 , 在 金融 科 ”基础 移动 方式 , 改善 了 算法 搜索 停滞 .陷入 局 部 最 优 的 缺陷 。 


EE 


Hr 


阔 的 发 展 前 景 。 在 使 用 反 向 学 习 框 架 改进 黏 菌 算法 的 研究 中 ， 其 大 多 使 用 单 
、 稳 定性 差 、 求 解 精 反 向 策略 ， 且 对 反 向 种 群 重 构 的 时 机 缺乏 判断 和 控制 ， 在 


算法 依旧 存在 着 
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录用 定稿 


保持 算法 原本 特征 和 减少 算法 时 间 复 杂 度 上 存在 进一步 优化 
空间 ; 在 改进 医 菌 算法 基本 移动 方式 的 研究 中 ， 其 主要 使 用 
经 典 改进 策略 和 其 他 种 类 算法 的 核心 移动 方式 奉 换 黏 菌 算法 
移动 方式 ， 对 秋 菌 算法 不 同 搜索 行为 的 阶段 缺乏 区 分 和 针对 
性 调整 ， 在 保持 算法 原 有 优势 的 基础 上 弥补 算法 局 限 性 上 存 
在 进一步 优化 空间 。 
针对 黏 菌 算法 搜索 停滞 、 稳 定性 差 的 问题 ， 本 文 使 用 反 
向 学 习 的 策 框架 改进 糙 苗 算法 并 在 传统 反 向 学 习 的 基础 
上 针对 黏 菌 算法 自身 特点 进行 策略 的 调整 和 组 合 ， 提 出 了 无 
迹 西 格 玛 点 引导 的 拟 反 向 黏 菌 算法 (UQSMA)。 使 用 拟 反 向 和 
拟 反 射 的 综合 反 向 过 程 跳出 局 部 最 优 ， 保 持 种 群 多 样 性 ;， 同 
时 根据 种 群 多 样 性 适时 使 用 反 向 种 群 重 构 原 始 种 群 ， 减 少 固 
定 反 向 过 程 的 运算 开销 ， 加 强 算法 的 稳定 性 ， 使 用 无 迹 西 格 
玛 点 引导 种 群 向 理论 最 优 逼 近 ， 针 对 黏 菌 算法 不 同 的 搜索 阶 
段 使 用 不 同 的 改进 策略 ， 弥 补 算法 在 面 对 不 同 实验 和 不 同 问 
题 时 稳定 性 较 差 的 局 限 。 多 种 策略 平衡 了 种 群 的 勘探 和 开 
采 阶 段 ， 各 策略 组 合 后 在 保证 单一 策略 有 效 的 基础 上 弥补 
了 单一 策略 的 局 限 。 实 验 中 ， 测 试 了 CEC2017 基准 测试 函 
数 ， 实 验 结果 验证 了 改进 策略 对 算法 求解 精度 和 和 鲁 棒 性 的 
积极 影响 。 
1 ， 基 本 黏 菌 算 法 

秋 菌 算法 中 的 黏 菌 一 般 是 指 多 头 绒 泡 菌 ， 因 其 最 初 被 归 
类 为 真菌 ， 所 以 被 命名 为 " 黏 菌 ”。 黏 菌 是 一 种 生活 在 寒冷 潮 
湿地 方 的 真 核 生 物 ， 寿 菌 的 活跃 动态 阶段 是 医 菌 算法 研究 的 
主要 阶段 。 在 这 个 阶段 ， 黏 菌 使 用 体内 的 有 机 物 寻找 食物 ， 


| 上 中 


包围 它 ， 并 分 泌 酶 来 消化 它 。 在 黏 菌 迁移 过 程 中 ， 其 前 端 会 
延伸 成 扇形 ， 后 端 则 连接 着 一 个 静脉 网 络 ， 保 证 旨 
中 流动 ， 从 而 达到 信息 的 沟通 ， 并 且 根 据 搜寻 到 食物 的 质量 
和 密度 动态 调整 静脉 直径 的 大 小 ， 最 终 将 整个 黏 菌 的 位 置 调 
整 至 最 优 食 物 处 。 受 到 黏 菌 这 种 行为 模式 的 启发 ，Shimin Li 
及 其 团队 设计 开发 了 黏 菌 算法 。 
基本 黏 菌 算法 的 抽象 数学 描述 如 下 : 假设 个 黏 菌 分 布 
在 4 维 搜索 空间 中 ,第 1 次 迭代 第 i 个 条 菌 的 位 置 可 以 表示 为 
Xi(1) ，max_f 为 最 大 迭代 次 数 , 全 局 最 优 位 置 X, 表示 食物 位 
置 。 黏 菌 的 砚 食 寻 优 过 程 分 为 三 个 阶段 : 
1) 接近 食物 
黏 菌 可 以 根据 空气 
X, (t+vb-(W: X(t)— X(t)),r<p 


x(e+)-{ ve:X(t),r>p (1) 
其 中 ， w 是 服从 [-a.4] 均匀 分 布 的 系数 ，vc 从 0 到 1 线性 递 
增 ，: 代表 当前 迭代 轮 次 ，X; 代表 目前 为 止 找 到 食物 浓度 最 
高 的 个 体 的 位 置 ， 即 当前 全 局 最 优 解 ，X 代表 黏 菌 位 置 ，X。 


中 的 气味 接近 食物 ， 下 式 模 拟 了 其 接 
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l+r: os( 好 1] condition 
W (SmellIndex(i))= 学 本 (5) 
1—r: os( + 小 others 
bF —wFr 
SmellIndex = sort(S) (6) 


其 中 选择 分 支 的 判断 条 件 condition 为 适应 度 排名 位 于 前 一 半 
种 群 ， 7 为 [0 间 的 一 个 随机 数 ， max :表示 最 大 友 代 次 数 ， 
bF 表示 当前 轮 次 获得 的 最 优 适应 度 值 ， 相 对 的 ， wF 表示 当 
前 轮 次 获得 的 最 差 适应 度 值 ， SmellIndex 表示 按照 适应 度 值 大 
小 排序 的 顾 菌 个 体 索 引 (在 求解 最 小 值 问 题 时 使 用 升序 排序 
方法 )。 

2) 缠绕 食物 

该 阶段 模拟 了 条 菌 搜索 食物 过 程 中 静脉 结构 的 收缩 模式 ， 
当 黏 菌 静 肪 接触 到 的 食物 浓度 越 大 ， 其 通过 生物 振荡 反应 产 
生 的 波 就 越 强 ， 细 胞 质 流动 就 越 快 ， 静 肪 细胞质 粘 稠 程度 也 
就 越 高 。 式 (53) 描 述 了 黏 菌 在 搜索 区 域 基于 食物 浓度 或 质量 的 
判断 获得 的 反馈 ， 若 该 位 置 的 食物 浓度 处 于 整个 种 群 截至 
前 所 有 搜索 结果 的 上 位 ， 则 给 予 该 位 置 务 菌 正 向 反馈 ， 反 之 
给 予 下 位 欠 菌 反 向 反馈 ， 以 此 方式 充分 利用 适应 度 的 计算 结 
果 信 息 来 调整 不 同位 置 黏 菌 的 权重 ， 进 一 步调 整整 个 种 群 的 
搜索 方式 。 同 时 ， 随 机 数 " 保证 了 黏 菌 个 体 间 的 多 样 性 。 


基于 上 述 理论 ， 翻 菌 的 位 置 更 新 公式 如 下 ; 
rand (UB-—L1B)+LB,rand <z 
X*=4X,(t)+vb:(W: X(t)— Xs(t)),r<p (7) 
ve:X(1),r2p 


其 中 LB 和 UB 代表 了 搜索 空间 的 上 下 边界 约束 ，wand 和 7 代 
表 了 [0 间 随 机 数 。 转 换 概率 z 在 Shimin Li 的 参数 实验 中 被 
证 设置 为 0.03 时 整个 算法 表现 最 好 ， 故 本 文中 z 同样 设置 为 
固定 值 0.03 。 
3) 振荡 反应 
黏 菌 依 靠 生物 的 震荡 反映 来 调整 静脉 中 细胞 质 的 流动 从 
而 调整 自身 位 置 搜索 最 优 食物 位 置 。 黏 菌 算法 中 的 参数 允 、 
ww 、vc 用 来 模拟 秋 菌 的 振荡 反应 ，W 根据 适应 度 选择 不 同 振 
荡 模 式 的 行为 模拟 了 黏 菌 个 体 根据 其 他 位 置 个 体 搜索 到 的 食 
物 情况 选择 自身 下 一 步行 为 的 搜索 模式 ，w 动态 调整 振荡 过 
程 中 使 用 其 他 个 体 信 息 的 口径 模拟 了 黏 菌 静脉 直径 变化 的 过 
程 ，vc 线性 变化 模拟 了 黏 菌 对 自身 历史 信息 保留 程度 的 变化 。 


2 ”无 迹 西 格 玛 点 引导 的 拟 反 向 忒 菌 算法 


2.1 反 向 学 习 策略 

Tizhoosh0 首次 提出 了 反 向 学 习 (opposition-based 
learning, OBL) 的 概念 , 其 广泛 应 用 于 群 智能 算法 的 改进 策略 
中 ， 同 时 反 向 学 习 也 具有 许多 扩展 和 变 体形 态 。 反 向 学 习 策 
略 有 助 于 扩大 搜索 空间 ， 加 速算 法 收敛 ， 维 持 种 群 多 样 性 。 


和 XX， 代表 两 个 从 锋 菌 种 群 中 随机 抽取 的 个 体 ，w 代表 笑 苗 
个 体 的 权重 。 
计算 7 的 公式 如 下 : 
p=tanh|s (i)-DF| (2) 
其 中 iel2,.…,n ，5( 引 表示 黏 菌 个 体 的 适应 度 值 ，DF 表示 所 有 
迭代 轮 次 中 找到 的 最 优 适应 度 。 
计算 w 的 公式 如 下 : 
vb=[-a,a] (3) 
其 中 a 为 随和 迭代 次 数 减 小 的 动态 参数 ， 用 来 收缩 随机 数 w 的 
生成 范围 ， 以 此 模拟 黏 菌 接近 食物 源 过 程 中 静脉 逐渐 收缩 提 
高 寻 优 精度 的 过 程 ，2 的 计算 公式 如 下 : 


二 4) 


a=arctanh [- 


计算 权重 W 的 公式 如 下 : 


反 疝 学 习 的 核心 机 制 是 同时 评估 原始 种 群 和 反 向 种 群 的 适应 
度 ， 并 从 中 挑选 适应 度 更 优 的 个 体 组 成 新 的 种 群 参与 后 续 计 
算 。 传 统 反 向 点 生成 的 数学 表达 式 如 下 : 
x9 =LB; +UB,—xX; (8) 
拟 反 向 学 习 04(quasi-opposition-based learning, QOBL) 和 
拟 反 射 学 习 05](quasi-reflection-based learning，QRBL) 是 反 向 
学 习 的 两 种 变 体 ， 拟 反 向 和 拟 反射 点 生成 的 数学 表达 如 下 : 


LB,; +UB, 
2 -dl 3 J +UB, 的 (9) 
de | 
XY = ran 二 ,Xj (10) 


文献 [16] 使 用 了 哈里 斯 诬 算 法 (HHO) 中 控制 种 群 实施 不 
同 搜索 行为 的 参数 控制 通过 拟 反 向 和 拟 反 射 过 程 生 成 的 反 向 
种 群 。 证 明了 拟 反 向 学 习 过 程 可 以 规避 传统 反 向 学 习 位 置 转 
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换 过 大 导致 的 精度 不 高 问题 ， 同 时 可 以 相对 地 扩大 搜索 空间 
加 速 种 群 收敛 ， 拟 反射 学 习 可 以 扩大 个 体位 置 附近 的 搜索 空 
间 ， 提 高 种 群 的 搜索 精度 。 图 1 展示 了 一 维 空间 内 拟 反 向 学 
习 和 拟 反射 学 习 相 对 于 传统 反 向 学 习 的 搜索 空间 。 


ib 并 lbtub Vo ub 


图 1 三 种 反 向 学 习 在 一 维 空 间 的 搜索 空间 示意 图 
Fig.1 Search space for three OBL in 1-dimensional space 


标准 的 群 智能 算 法 框架 中 存在 着 明确 的 勘探 和 开采 行为 


从 而 在 种 群 
干扰 。 
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陷入 停滞 时 进行 帮助 ， 在 种 群 快速 收敛 时 不 进行 


2.3 ”无 迹 西格玛 点 引导 策略 
文献 [20] 首 次 提出 了 使 用 无 迹 变 换 (unscented 
transformation，UT) 生 成 一 系列 无 迹 西格玛 点 帮助 原始 种 群 


晕 近 型 


E 论 最 优点 位 置 。 通 过 修正 原始 算法 的 步 进 方式 ， 丰 富 


种 群 搜 索 的 方式 ， 提 高 算法 面 对 不 同 问题 时 的 泛 用 性 ， 避 免 


算法 陷 


入 停滞 ， 无 迹 西格玛 点 数学 表达 式 如 下 : 


G0 = 工 
=Xs+0g4,d =1,2,...,D (14) 
Gaip =X 0,d =1,2,...,D 


的 区 分 ， 若 使 用 单一 的 反 向 学 习 策 略 ， 生 成 的 反 向 种 群 特征 


中 忆 为 维度 4 上 位 置 的 平均 值 ,0 为 维度 4 上 位 置 的 标准 


由 较为 单一 ， 对 算法 勘探 和 开采 阶段 的 辅助 效果 会 不 同 ， 且 
会 存在 一 个 效果 相对 不 好 的 阶段 。 为 了 更 有 效 地 构建 反 向 种 
群 ， 本 文 在 黏 菌 算法 进行 勘探 时 使 用 大 步 长 粗糙 搜索 的 拟 反 
名 学习 过 程 生成 反 向 个 体 ， 在 开采 时 使 用 小 步 长 精细 搜索 的 


使 用 


由 上 式 可 知 无 迹 西格玛 点 共 2D+1 个 。 本 文 使 用 第 一 个 无 
迹 西格玛 点 , 即 种 群 位 置 的 均值 来 引导 秋 菌 算法 的 勘探 阶段 ; 
剩余 无 迹 西格玛 点 ， 即 使 用 标准 差 和 均值 运算 得 出 的 新 


拟 反射 学 习 过 程 生成 反 向 个 体 ， 最 终生 成 完整 的 综合 反 向 种 
群 ， 提 升 反 向 学 习 框架 针对 不 同 搜索 空间 特征 的 优化 效果 。 
在 黏 菌 算法 中 ， 参 数 控制 着 种 群 的 勤 探 和 开采 行为 ， 本 文 
在 医 菌 进行 勘探 时 使 用 拟 反 向 过 程 生 成 反 向 个 体 ， 在 年 医 菌 
进行 开采 时 使 用 拟 反 射 过 程 生成 反 向 个 体 ,数学 表达 式 如 下 : 


rand (HE (8+U8-n) |r< p 
OX = (11) 


rand (2 ),x lr> p 


其 中 参数 Pp 随 着 从 代 轮 次 的 增加 而 减 小 ， 即 生成 的 综合 反 向 
种 群 在 算法 前 期 以 拟 反 向 个 体 为 主 ， 在 算法 后 期 以 拟 反 射 个 
体 为 主 ， 从 而 动态 平衡 反 向 种 群 内 个 体 勘 探 和 开采 的 数量 比 
例 。 
2.2 线性 多 样 性 度量 


数学 表达 式 如 下 
D, (2)=- BY (a (12) 


CEF- 训 2ot7 到- 让 2 
) = 区 三 zf 13 
NP 守 NP 和 (13) 


其 中 ，D 为 问题 维度 ，4 为 当前 维度 ，WNP 为 种 群 规模 。 
传统 反 向 学 习 中 ， 每 轮 迭 代 生 成 的 反 向 种 群 都 会 与 原始 
种 群 一 起 进行 比较 从 而 生成 新 种 群 ， 不 仅 会 大 大 提升 运算 时 
间 开 销 ， 新 种 群 位 置 的 重 置 还 会 打 断 原始 种 群 的 收敛 过 程 ， 
导致 算法 在 多 次 独立 实验 时 表现 出 不 稳定 的 问题 ， 因 此 适时 
地 进行 反 向 重 构 种 群 才 能 够 最 大 程度 的 发 挥 反 向 学 习 策略 的 
优势 。 

受 文献 [19] 启 发 ， 本 文 设 置 一 个 种 群 多 样 性 阔 值 
(diversity threshold, DT) 用 以 调控 当前 种 群 是 否 需 要 通过 上 述 
综合 反 向 种 群生 成 新 种 群 。 若 计算 得 出 的 种 群 多 样 性 大 于 阐 
值 ， 则 生成 新 种 群 ， 若 种 群 多 样 性 小 于 立 值 ， 则 不 生成 新 种 
群 。 由 于 笑 菌 算法 本 身 具 有 维持 自身 种 群 多样 性 的 能 力 ， 故 
当 种 群 多样 性 变 小 时 可 能 表明 种 群 处 于 收敛 过 程 。 在 实验 中 
发 现 ， 即 使 颖 菌 算法 陷入 局 部 最 优 ， 其 种 群 多 样 性 并 不 会 始 
终 维 持 低 水 平 ， 而 是 会 在 高 低 水 平 间 震 荡 波 动 ， 此 时 车 种 群 
多 样 性 超过 多 样 性 闵 值 ， 综 合 反 向 种 群 便 可 以 通过 生成 新 种 
群 帮助 原始 种 群 扩 大 搜索 范围 ， 跳 出 局 部 最 优 。 该 机 制 以 较 
小 的 计算 量 达 到 了 缩减 反 向 重 构 过 程 的 效果 ， 显 著 的 减少 了 


文献 [17] 提 出 了 度量 进化 算法 中 种 群 多 样 性 的 度量 方式 ， 


位 置 来 引导 恬 菌 算法 的 开采 阶段 ， 帮 助 算法 在 反 向 学 习 扩展 
后 的 搜索 空间 中 


核 ， 


| 导 算法 逼近 理论 最 优 。 式 (7) 为 香 菌 算法 的 


心 搜索 框架 ， 其 中 <P 时 种 群 表现 出 勘探 的 行为 ， 其 基本 


移动 方式 为 使 用 最 优 位 置 部 分 信息 的 基础 上 使 用 两 个 随机 个 
体位 置信 息 生 成 步 长 , 本 文 将 X 蔡 换 为 X， 减少 步 长 生成 的 


随机 性 ， 使 用 均值 引导 种 群 在 勘探 时 快速 收敛 ;，r>p 时 种 群 


表现 出 开采 的 行为 ， 其 基本 移动 方式 为 对 自身 位 置 的 衰减 ， 


本 文 将 其 蔡 换 为 均值 与 方差 生成 的 新 位 置 , 同时 保留 参数 ve 
保证 一 定 的 随机 性 ， 帮 助 斐 菌 算法 开采 阶 段 具 有 荧 近 理论 


最 优点 的 能 


， 同时 设 定 转换 概率 , 不 全 部 计算 2D 个 无 迹 西 


格 玛 点 ， 节 省 运算 开销 。 本 文 使 用 无 迹 西 格 玛 点 改进 后 的 式 
(7) 如 下 : 


rand (UB-LB)+LB,rand <z 
X(t)+vb:(W-X() -X(t)),r<p 
X(t) +ve:o(t),r>p r'>0.5 
X(t) veo(t),r2p r'<0.5 


(15) 


2.4 无 迹 西 格 玛 点 引导 的 拟 反 向 黏 菌 算法 


本 文 提 出 了 


无 迹 西格玛 点 引导 的 拟 反 向 甜 菌 算法 


(unscented sigma point guided quasi-opposite Slime Mould 
Algorithm，UQSMA)。 首 先 ， 根 据 原始 条 菌 算法 勘探 和 开采 


的 不 同 阶段 ， 使 


j 拟 反 向 学 习 和 拟 反 射 学 习 过 程 生成 综合 反 


名 种 群 ， 在 番 菌 种 群 不 同 的 空间 特征 搜索 阶段 使 用 不 同 的 反 


句 过 程 有 助 于 放大 黏 苗 在 不 同 搜索 行为 时 的 倾向 ， 强 化 黏 攻 


算法 本 身 的 行为 特征 ， 加 速算 法 的 启动 ， 其 次 ， 使 用 线性 多 
样 性 度量 计算 种 群 多 样 性 ， 并 根据 种 群 的 实时 多 样 性 判断 是 


度 ， 


格 玛 点 蔡 换 原始 笑 


否 使 用 综合 反 向 利 


h 群 加 入 原始 种 群生 成 新 种 群 ， 保 证 算法 精 


节约 运算 开销 ， 提 高 算法 的 稳定 性 ， 最 后 ， 使 用 无 迹 西 


菌 算 法 的 基本 移动 方式 ， 引 导 种 群 逼 近 理 


论 最 优 解 ， 丰 富 搜索 的 路 径 ， 使 算法 保持 最 优 解 的 更 新 和 优 


化 ， 避 免 算 法 停 光 。 


2.5 改进 算法 的 步骤 法 

UQSMA 算法 的 流程 图 如 图 2 所 示 。 

2.6 ”时 间 复 杂 度 分 析 

假设 种 群 规模 为 NP , 搜索 空间 维度 为 D ,最 大 和 迭代 次 数 


为 7 ， 原 始 SMA 算法 的 时 间 复 杂 度 


E 要 来 源 为 初始 化 、 适 


应 度 计算 和 位 置 更 新 ， 其 时 间 复 杂 度 为 
O(SMA)=O(NP)+O(NP*T)+O(NP*D*T) (16) 
UQSMA 算法 主要 由 拟 反 向 和 拟 反 射 学 习 组 成 的 综合 反 


句 学习 


运算 的 时 间 开销 ;同时 由 于 种 群 多 样 性 与 种 群 的 收敛 与 否 
有 较 强 相关 性 ， 因 此 通过 设 定 阔 值 对 种 群 状 态 进 行 筛 选 后 习 


/4 


几 


构 得 到 的 新 种 群 也 一 定 程度 上 可 以 对 应 种 群 的 不 同 收敛 状态 ， 


喇 


NP 个 反 向 个 体 的 
成 新 种 群 , 该 过 程 


匡 架 和 无 迹 西格玛 点 引导 机 制 两 个 改进 策略 组 成 。 反 
过 程 需要 根据 NP 个 个 体 生 成 NP 个 反 向 个 体 ， 增 加 计算 


适应 度 ， 对 所 有 2NP 个 个 体 进 行 排序 并 生 


时 间 复 杂 度 为 O(NP*D+2NP+2NPlg2NP), 与 
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标准 反 向 学 习 算法 相同 ; 无 迹 西格玛 点 生成 过 程 不 增加 时 间 定义 为 

复杂 度 ， 则 UQSMA 算法 的 时 间 复 杂 度 如 下 : DGx) | Gr f(s(X)>7(G) 
O(UQSMA)=O(NP+(NP*D+2NP+2NPlg2NP)) ~ lg(X,), f(g(Xi))<f(G.) 
+O(NP*T)+O((2NP*D+2NP+2NPlg2NP)*T) (17) G， f(n(xX,))zf(G,) 
=O(NP*(3+ D+2lg2NP+T(3+2D+21g82NP))) 和 DP(G,X)= A(Xi), f(A(Xi)) <f(G,) 


UQSMA 算法 实际 运行 时 ， 由 于 存在 适应 度 病 值 判断 是 中 ，8 对 应 拟 反 向 与 拟 反 射 学 习 过 程 ， 与 标准 反 向 学 习 过 
否 进 行 综合 反 向 过 程 ， 故 根据 实际 问题 的 不 同 算法 总 时 间 复 。 程 在 位 置 更 新 步骤 上 逻辑 一 致 ， < 表示 R* 中 第 1 代 个 体 i 的 
杂 性 有 不 同 程度 的 降低 ， 即 0((2NP*D+2NP+2NPlg2NP)*T) 部 样本 空间 , 8 (Xi) 表 示 第 1 次 迭代 中 通过 拟 反 向 和 拟 反 射 学 习 
分 在 实际 运算 过 程 中 会 更 小 ， 故 式 (17) 是 UQSMA 算法 时 间 ”生成 的 综合 反 向 种 群 与 原始 黏 菌 种 群 整合 后 更 新 的 个 体 位 


并 


复杂 度 的 理论 最 大 值 。 置 ，A(X) 表示 黏 菌 个 体 : 经 过 位 置 移动 后 第 上 次 更 新 后 的 位 
置 ，G, 为 当前 得 到 的 全 局 最 优 位 置 。 由 上 述 表 达 式 可 知 ，C 
开始 按 式 (12) 计 算 所 对 应 的 全 局 最 优 适 应 度 单 调 不 增 ， 且 逐渐 向 解 空间 的 下 确 
种 群 线性 多 样 性 界 收敛 。 
证 明成 立 。 
Les 多 样 性 是 否 天 ~ 二 定理 2 UQSMA 满足 条 件 2。 
| Ed 证 明 车 算法 满足 条 件 2, 种 群 规模 为 NP 的 入 菌 种 群 的 
We 和 样本 空间 必须 包含 s ， 即 
计算 SMA 种 群 和 综 计算 SMA ge U i 
计算 SVMA 种 群 和 反 合 反 向 种 群 适应 度 | | 。 种 群 适应 度 Te 
| RH 其 中 ， Mi, 为 第 + 次 迭代 中 个 体 i 的 样本 空间 支撑 集 。 
< 不 肠 全 荆 ua 、 人 
生成 新 种 群 对 于 UQSMA 中 拟 反 向 和 拟 反射 学 习 过 程 的 位 置 更 新 机 
二 全 和 过 制 ， 有 
| 按 式 (15) 进 行 有 
位 置 更 新 rand| (EE (LB+UB -xX, )| < p 
策 西 全 局 最 优 位 置 8 LB+UB 
和 最 优 适应 度 按 式 (1) 生 成 md 人 /| 
综合 反 向 种 群 


n 


拟 反 向 和 拟 反 射 学 习 过 程 后 Xi; 到 达 s (Xi) 的 概率 大 于 零 。 
图 2 UQSMA 算法 流程 图 此 随 着 迭代 次 数 的 增加 每 个 8(Xi,) 的 闭 包 V[s(X)] 与 其 并 集 


> ff Fig.2 Flow chart ofUQSMA 六 up N | ,到 扩 3 
CN 27 收 全 性 分 析 Us (Xi,) 的 闭 包 V[( 放 8(X,)| 都 在 逐渐 扩大 ， 因此 存在 一 人 


UQSMA 算法 本 质 是 一 种 随机 搜索 算法 ， 可 由 概率 测度 。 整数 4 ， 当 1>4 时 

法 证 明 UQSMA 的 全 局 收敛 性 。 全 局 收敛 性 准则 及 定理 ba 闻 本 

述 了 算法 满足 全 局 收敛 性 需要 满足 的 两 个 条 件 ; 

条 件 1 7(D(z6)s7()， 且 若 ses， 则 辣 理 ， 原 始 务 菌 算法 的 位 置 更 新 过 程 得 到 的 样本 空间 并 

ee 集 的 闭 包 Y[J"n(x.)] 也 在 逐渐 扩大 。 因 此 存在 一 个 整数 

中 ， 为 适应 度 函 数 ， 为 随机 搜索 算法 ，: 为 解 空间 以 Ce 

子 集 $ 的 其 中 一 个 点 ， 该 点 能 使 适应 度 函数 取得 最 小 值 或 在 。 当 1>+ 时 

s 上 产生 可 接受 的 函数 值 下 确 界 ，& 为 算法 D 在 交代 搜索 过 

程 中 得 到 的 解 。 
条 件 2 对 于 s 中 任意 的 Borel 子 集 4， 车 其 概率 测度 

V[4]>0 3 则 有 


此 可 知 8(X) 为 R* 中 不 同 于 的 样本 空间 ， 即 经 过 
天 
集 


5 


U 


(Ur) Us) 


令 s 的 Borel 子 集 D=(UA(Xi))J(Uis(X)) ， 有 


IIa-w(4)-0 VID]>0 ，A(D)=Y"A(D)=1， 即 


其 中 ，44 (4)=P(x EA,?,…x 中 ) 为 算法 D 第 k 次 迭代 结果 在 
眼 合 4 上 的 概率 测度 ， {x%}。 该 假设 表明 ， 对 于 5 中 任意 满 


[1-4tpD=0 


足 Y[4]>0 的 子 集 4 经 过 算法 的 无 穷 次 迭代 后 , 错过 解 空间 5 证 明成 立 。 

的 任意 的 Borel 子 集 4 的 概率 为 0， 即 满足 上 述 条 件 的 算法 综 上 ， 由 引 理 1 得 ，UQSMA 具有 全 局 收敛 性 ， 且 由 于 

在 无 穷 次 搜索 后 无 法 达到 近似 全 局 最 优点 的 概率 为 0。 拟 反 向 和 拟 反 射 学 习 过 程 生 成 的 新 样本 空间 不 同 于 原始 慕 菌 
引 理 1 随机 优化 算法 全 局 收 和 敛 定理 08. 若 函数 /可 测 ， ”种 群 的 样本 空间 ， 该 策略 可 以 扩大 整个 种 群 在 一 次 迭代 中 搜 


搜索 空间 5 是 R* 上 的 可 测 子 集 , 若 算 法 满足 条 件 1 和 条 件 2， 索 的 样本 空间 , 故 收敛 速度 相对 于 原始 黏 菌 算法 定 加 快 。 
则 有 本 
limP[xe eR]=1 3 ”仿真 实验 及 结果 分 析 

其 中 ，R 是 全 局 最 优点 集合 ，P[* eR] 是 算法 第 上 次 搜索 欠 代 实验 环境 为 Windows10，64 位 操作 系统 ，CPU 为 AMD 
结果 落 在 R 中 的 概率 。 Ryzen 7 5800H， 主 频 3.20GHz， 内 存 16GB ， 算 法 基于 

定理 1 UQSMA 满足 条 件 1。 MATLAB2020b 编写 。 

证 明 ”本 文 算法 最 优 位置 更 新 包含 原始 黏 菌 算法 的 位 置 实验 一 “为 了 充分 验证 本 文 所 提出 的 UQSMA 在 求解 函 
更 新 和 拟 反 向 与 拟 反 射 学 习 的 位 置 更 新 ， 故 可 将 算法 D 抽象 。” 数 优化 问题 时 的 求解 精度 和 稳定 性 ， 将 UQSMA 算法 与 其 他 
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本 文 改进 策略 能 够 使 算法 保证 持续 的 搜索 效果 和 搜索 效率 ， 
克服 原始 黏 菌 算法 搜索 停止 的 局 限 。 在 此 基础 上 UQSMA 算 
法 的 求解 精度 均 高 于 传统 反 向 改进 的 黏 菌 算法 算法 
DOLSMA, 说 明 拟 反 向 和 拟 反 射 学 习 的 组 合 对 于 顾 菌 算法 效 
果 的 提升 优 于 传统 的 动态 反 向 学 习 。 同 时 ，UQSMA 算法 的 
方差 整体 低 于 其 他 算法 ， 说 明了 本 文 策略 的 组 合 对 于 算法 稳 
定性 的 提升 是 有 效 的 ， 对 于 算法 最 终 求解 精度 和 重 棒 性 具有 
一 定 的 效果 。 


T 


算法 在 CEC2017 测试 函数 上 进行 求 进行 实验 结果 的 对 
比 。 对 比 用 算法 包括 最 新 的 改进 颖 菌 算法 AFSMAD1、 动 态 反 
向 改进 的 黏 菌 算 法 DOLSMA、 精 英 反 向 与 其 他 策略 结合 改进 
的 鲸鱼 优化 算法 EGSWOARPR34、 原 始 的 寿 菌 算法 SMAII、 同 
科研 团队 提出 的 哈里 斯 鹰 算 法 HHOP3 和 经 典 的 鲸鱼 优化 
算法 WOAPH。 
为 保证 每 种 算法 对 比 时 的 公平 性 ， 实 验 中 种 群 规 模 设 置 
为 30， 最 大 迭代 次 数 设置 为 500, 并 独立 运行 30 次 计算 得 出 
统计 结果 。 对 比 用 算法 参数 设置 均 来 源 于 对 应 文献 设置 ， 如 
表 1 所 示 。 


不 | 


人 


表 1 
Tab. 1 


主要 参数 


Main parameters 


3) 


本 文选 用 了 表 2 所 列 示 的 29 个 CEC2017 基准 测试 函数 ， 算法 主要 参数 


其 类 型 包含 单 
合 (Hybrid, H) 及 复合 (Composition, C)， 旨 
50、100， 此 处 实验 使 用 标准 旨 
的 不 稳定 性 习惯 上 不 使 用 。 

表 3 为 UQSMA 与 第 
总 体 上 看 ，UQSMA 在 29 个 CEC2017 基准 测 
个 函数 上 平均 求解 精度 的 表现 优 于 其 他 所 有 


茜 (Unimodal, UM)、 多 峰 (Multimodal, MMD)、 混 
住 度 可 选用 10、30、 
度 30 维 ，02 号 函数 由 于 高 维 


UQSMA 
AFSMA 
DOLSMA 
EGSWOA 
SMA 
HHO 
WOA 


z=0.03 ; DT=0.3 
o=1 
Z=0.03 
b=1 ; x1=-nx+(1-7)*2n ; x2=-n+t*27n 
Z=0.03 
B=1.5 
b=1 


部 分 算法 进行 对 比 的 实验 结果 。 
试 函数 中 在 21 
对 比 算法 ， 说 明 


表 2 CEC2017 函数 
Tab.2 CEC2017 function 


函数 函数 名 称 类 型 最 优 值 | 函数 函数 名 称 类 型 最 优 值 


CEC01 
CEC03 
CEC04 
CEC05 
CEC06 
CEC07 
CEC08 
CEC09 
CEC10 
CEC11 
CEC12 
CEC13 
CEC14 
CEC15 


100 
300 
400 
500 
600 
700 
800 
900 
1000 
1100 
1200 
1300 


CEC16 
CEC17 
CEC18 
CEC19 
CEC20 
CEC21 
CEC22 
CEC23 
CEC24 
CEC25 
CEC26 
CEC27 


1600 
1700 
1800 
1900 
2000 
2100 
2200 
2300 
2400 
2500 
2600 
2700 
2800 
2900 


Shifted and Rotated Bent Cigar Function 
Shifted and Rotated Zakharov Function 
Shifted and Rotated Rosenbrock’s Function 
Shifted and Rotated Rastrigin’s Function 
Shifted and Rotated Expanded Scaffer’s F6 Function 
Shifted and Rotated Lunacek Bi-Rastrigin Function 
Shifted and Rotated Non-Continuous Rastrigin’s Function 
Shifted and Rotated Lévy Function 
Shifted and Rotated Schwefel’s Function 
Hybrid Function 1 (N=3) 

Hybrid Function 2 (N=3) 

Hybrid Function 3 (N=3) 


Hybrid Function 5 (N=4) 
Hybrid Function 6 (N=4) 
Hybrid Function 6 (N=5) 
Hybrid Function 6 (N=5) 
Hybrid Function 6 (N=6) 
Composition Function 1 (N=3) 
Composition Function 2 (N=3) 
Composition Function 3 (N=4) 
Composition Function 4 (N=4) 
Composition Function 5 (N=5) 
Composition Function 6 (N=5) 
Composition Function 7 (N=6) 
Hybrid Function 4 (N=4) 1400 |CEC28 Composition Function 8 (N=6) 
Hybrid Function 5 (N=4) 1500 |CEC29 Composition Function 9 (N=3) 


和 之 之 之 之 之 Zcc 
和 REERERRZE 
人 人 


表 3 CEC2017 测试 函数 和 Wilcoxon 秩 和 检验 的 结果 对 比 
Tab.3 Comparison of CEC2017 function results and Wilcoxon’s rank-sum test 


UQSMA AFSMA DOLSMA EGSWOA SMA HHO WOA 


Mean(Std)win 


Mean(Std)win 


Mean(Std)win 


Mean(Std)win 


Mean(Std)win 


Mean(Std)win 


Mean(Std)win 


CEC01 
CEC03 
CEC04 
CEC05 
CEC06 
CEC07 
CEC08 
CEC09 
CEC 
CEC! 
CEC 
CEC 
CE 
CEC 
CEC 
CEC 
CEC 
EC 
CEC20 
CEC21 
CEC22 
CEC23 
CEC24 
CEC25 
CEC26 
CEC27 
CEC28 
CEC29 
CEC30 
+/=/- 


oo 和 ww 一 口 


1.21E+04(8.85E+03) 
1.26E+04(5.16E+03) 
$5.01E+02(2.58E+01) 
6.04E+02(3.07E+01) 
6.19E+02(7.99E+00) 
8.S1E+02(2.78E+01) 
9.00E+02(1.96E+01) 
1.99E+03(6.54E+02) 
4.79E+03(6.36E+02) 
1.26E+03(4.50E+01) 
3.56E+06(2.84E+06) 
6.82E+04(3.64E+04) 
4.07E+04(3.79E+04) 
4.58E+04(2.43E+04) 
2.61E+03(2.28E+02) 
2.13E+03(1.75E+02) 
1.02E+06(7.54E+05) 
2.45E+04(2.25E+04) 
2.53E+03(1.56E+02) 
2.39E+03(4.08E+01) 
3.0SE+03(1.53E+03) 
2.76E+03(2.58E+01) 
2.94E+03(4.54E+01) 
2.90E+03(1.71E+01) 
4.3SE+03(9.63E+02) 
3.24E+03(1.95E+01) 
3.24E+03(2.49E+01) 
3.9SE+03(2.02E+02) 
4.28E+05(2.27E+05) 


4.33E+04(1.07E+05)= 
4.80E+04(1.44E+04)= 
5.11E+02(2.75E+01)+ 
6.74E+02(3.36E+01)= 
6.24E+02(1.02E+01)+ 
9.35E+02(5.40E+01)+ 
9.54E+02(3.33E+01)+ 
5.05E+03(1.56E+03)+ 
5.24E+03(7.21E+02)+ 
1.2SE+03(4.85E+01)- 
3.37E+06(3.13E+06)= 
3.78E+04(6.34E+04)= 
5.05E+05(4.93E+05)+ 
1.33E+04(1.33E+04)- 
2.64E+03(3.03E+02)+ 
2.45E+03(2.64E+02)= 
2.45E+06(2.72E+06)+ 
1.19E+04(1.56E+04)- 
2.71E+03(1.88E+02)+ 
2.45E+03(3.67E+01)+ 
6.13E+03(1.85E+03)+ 
2.80E+03(3.96E+01)+ 
2.98E+03(4.81E+01)+ 
2.91E+03(2.90E+01)+ 
5.37E+03(4.12E+02)= 
3.25E+03(2.09E+01)+ 
3.27E+03(2.65E+01)= 
4.02E+03(2.94E+02)+ 
3.28E+04(1.43E+04)= 
17/9/3 


5.50E+05(3.24E+05)+ 
3.27E+04(1.16E+04)= 
5.03E+02(1.54E+01)+ 
6.37E+02(3.44E+01)= 
6.21E+02(1.16E+01)+ 
8.93E+02(4.03E+01)= 
9.31E+02(3.10E+01)+ 
4.39E+03(1.60E+03)+ 
5.05E+03(8.10E+02)+ 
1.28E+03(5.74E+01)= 
5.58E+06(4.07E+06)= 
1.01E+05(5.07E+04)+ 
1.79E+05(1.92E+05)+ 
2.78E+04(1.98E+04)+ 
2.61E+03(3.02E+02)+ 
2.36E+03(2.37E+02)= 
2.28E+06(3.35E+06)+ 
2.14E+04(1.97E+04)= 
2.60E+03(1.76E+02)= 
2.43E+03(3.34E+01)= 
5.82E+03(1.74E+03)+ 
2.77E+03(2.60E+01)+ 
2.96E+03(3.77E+01)= 
2.90E+03(1.44E+01)= 
4.83E+03(4.63E+02)= 
3.24E+03(1.61E+01)= 
3.27E+03(3.93E+01)= 
4.06E+03(2.25E+02)+ 
2.56E+05(2.14E+05)= 
14/15/0 


2.90E+10(6.55E+09)+ 
8.39E+04(6.63E+03)= 
6.34E+03(1.50E+03)+ 
8.99E+02(4.09E+01)+ 
6.83E+02(1.21E+01)+ 
1.38E+03(6.75E+01)+ 
1.12E+03(2.63E+01)+ 
1.08E+04(1.32E+03)+ 
8.39E+03(6.40E+02)+ 
7.12E+03(1.88E+03)+ 
2.02E+09(8.57E+08)+ 
4.07E+08(3.97E+08)+ 
2.36E+06(1.61E+06)+ 
7.88E+07(8.33E+07)+ 
4.62E+03(5.33E+02)+ 
2.92E+03(3.67E+02)+ 
1.84E+07(1.82E+07)+ 
9.63E+07(9.51E+07)+ 
2.90E+03(1.85E+02)+ 
2.67E+03(5.12E+01)+ 
7.88E+03(1.26E+03)+ 
3.24E+03(1.30E+02)+ 
3.44E+03(1.16E+02)+ 
3.85E+03(1.72E+02)+ 
9.24E+03(7.12E+02)+ 
3.20E+03(1.80E-09)- 
3.30E+03(4.00E-10)+ 
5.87E+03(7.31E+02)+ 
2.11E+08(1.30E+08)+ 
27/1/1 


1.19E+05(4.92E+04)+ 
3.61E+04(1.60E+04)= 
5.20E+02(3.98E+01)+ 
6.29E+02(3.93E+01)= 
6.18E+02(1.14E+01)- 
8.91E+02(4.95E+01)= 
9.34E+02(3.67E+01)+ 
4.37E+03(1.47E+03)+ 
4.97E+03(6.99E+02)+ 
1.29E+03(6.31E+01)= 
4.74E+06(3.82E+06)= 
6.60E+04(3.35E+04)= 
2.10E+05(2.20E+05)= 
2.47E+04(1.64E+04)+ 
2.56E+03(3.47E+02)- 
2.38E+03(2.18E+02)= 
2.20E+06(2.35E+06)+ 
1.77E+04(1.91E+04)+ 
2.65E+03(1.98E+02)= 
2.44E+03(3.85E+01)+ 
6.26E+03(9.23E+02)+ 
2.78E+03(3.14E+01)+ 
2.96E+03(3.56E+01)+ 
2.91E+03(2.34E+01)= 
5.02E+03(5.90E+02)= 
3.24E+03(1.77E+01)+ 
3.27E+03(4.50E+01)= 
4.04E+03(2.38E+02)+ 
1.17E+05(9.64E+04)= 
14/13/2 


3.81E+08(2.39E+08)+ 
5.88E+04(6.44E+03)= 
7.20E+02(9.25E+01)+ 
7.69E+02(3.14E+01)+ 
6.68E+02(5.38E+00)+ 
.32E+03(6.42E+01)+ 
9.93E+02(2.96E+01)+ 
8.92E+03(8.65E+02)+ 
6.09E+03(8.64E+02)+ 
.54E+03(1.70E+02)+ 
9.00E+07(1.01E+08)+ 
9.18E+06(3.98E+07)+ 
.40E+06(1.62E+06)+ 
.36E+05(9.08E+04)+ 
3.70E+03(5.09E+02)+ 
2.71E+03(3.45E+02)+ 
2.40E+06(1.92E+06)+ 
2.09E+06(1.70E+06)+ 
2.81E+03(1.87E+02)+ 
2.60E+03(5.07E+01)+ 
7.57E+03(1.01E+03)+ 
3.32E+03(1.35E+02)+ 
3.50E+03(1.80E+02)+ 
3.02E+03(3.77E+01)+ 
8.17E+03(8.36E+02)+ 
3.69E+03(1.90E+02)+ 
3.49E+03(1.24E+02)+ 
4.97E+03(5.31E+02)+ 
2.46E+07(4.10E+07)+ 
28/1/0 


5.65E+09(2.26E+09)+ 
2.80E+05(6.03E+04)= 
1.38E+03(3.85E+02)+ 
8.59E+02(5.33E+01)+ 
6.84E+02(1.31E+01)+ 
1.32E+03(7.53E+01)+ 
1.07E+03(4.85E+01)+ 
1.06E+04(3.37E+03)+ 
7.37E+03(6.18E+02)+ 
1.02E+04(4.98E+03)+ 
4.58E+08(2.46E+08)+ 
1.47E+07(2.11E+07)+ 
3.11E+06(2.72E+06)+ 
9.31E+06(1.94E+07)+ 
4.28E+03(5.92E+02)+ 
2.91E+03(2.84E+02)+ 
1.33E+07(1.27E+07)+ 
1.58E+07(1.31E+07)+ 
2.92E+03(2.15E+02)+ 
2.66E+03(5.08E+01)+ 
7.82E+03(1.84E+03)+ 
3.16E+03(9.95E+01)+ 
3.30E+03(9.78E+01)+ 
3.21E+03(8.34E+01)+ 
8.59E+03(1.14E+03)+ 
3.51E+03(1.51E+02)+ 
3.97E+03(4.12E+02)+ 
5.34E+03(5.61E+02)+ 
7.72E+07(5.38E+07)+ 
28/1/0 


E 改 进 算法 的 有 效 拆 


E， 需 要 对 群 智能 算法 
结果 进行 统计 检验 P4， 本 文选 择 Wilcoxon 秩 和 检验 进行 统 
生 水 平 设置 为 5%， 表 3 展示 了 UQSMA 算法 
每 列 最 后 一 个 符号 为 UQSMA 算法 


他 算法 的 检验 结果 ， 


| Wilcoxon 秩 和 检验 订 


异 , 即 UQSMA 算法 综合 优 于 对 比 算法 ; ”代表 无 论 均值 优 


值 结 果 优 于 


劣 且 Wilcoxon 秩 和 检验 证 明 两 算法 结果 分 布 不 


UQSMA 算法 综合 与 对 比 算法 相似 ;“- 代表 均值 结果 
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文 策略 间 组 合 可 以 达到 在 不 改变 原始 SMA 算法 搜索 特征 的 
基础 上 帮助 算法 持续 
算法 求解 精度 ， 
CEC2017 基准 


搜索 ， 防 止 算法 停滞 ， 从 而 提高 最 终 的 


可 以 证 明 UQSMA 的 持续 搜索 效果 在 求解 
则 试 函数 上 表现 更 优 。 

为 了 形象 的 观察 改进 算法 的 收敛 情况 与 稳定 性 ， 绘 制 
E 明 两 算法 结果 分 布 具有 差 ”” UQSMA 与 上 述 对 比 算法 的 收敛 曲线 和 箱 线 图 。 图 3 给 出 部 
分 基准 测试 函数 的 收敛 曲线 图 ,各 算法 曲线 图 例 与 图 3(h) 一 


有 明显 差 。 致 , 为 了 便于 观察 高 精度 的 曲线 收敛 情况 , 纵 坐 标 取 以 10 为 


劣 于 比较 算法 上 且 Wilcoxon 秩 逢 
有 差异 ， 即 UQSMA 算法 综合 劣 了 


E 明 两 算法 结果 分 布 
F 对 比 算法 。 最 后 一 行为 优 
合计 。 表 3 结果 表明 UQSMA 算 
结果 完全 不 劣 于 
相 较 于 改进 算法 
数 处 于 劣势 和 在 1 个 测试 函数 上 搜索 停滞 。 
表现 相似 ， 其 他 函数 均 处 于 优势 ， 相 较 于 AFSMA 算法 和 原 
始 SMA 算法 ，UQSMA 算法 仅 在 少数 测试 函数 上 处 于 劣势 。 


法 对 于 每 个 CEC2017 基准 测试 函数 的 求 角 
DOLSMA 、HHO 、WOA 三 个 对 比 算法 ; 
EGSWOA 仅 在 1 个 测试 函 


底 的 对 数 绘 表 
均值 的 箱 线 图 
后 期 保持 持续 搜索 的 
在 前 中 期 扩大 了 搜索 空间 ， 从 而 使 算法 快速 的 找到 更 优 的 寻 


所 有 曲线 图 像 。 图 4 给 出 部 分 基准 测试 函数 带 
。 观察 可 得 ，UQSMA 算法 具有 前 期 收敛 快 、 


优势 ， 说 明了 拟 反 向 过 程 和 拟 反射 过 程 


此 基础 上 最 终 的 求解 结果 多 数 优 于 原始 SMA 算法 ， 说 明 本 


优 方 向 从 而 力 
正 帮 助 算法 在 迭代 后 


分 析 可 知 UQSMA 算法 通过 设 定 种 群 多 样 性 阔 值 适时 的 进行 ”根据 多 样 性 进行 


; 无 迹 西格玛 点 对 于 种 群 搜索 行为 的 修 


期 依旧 能 够 保持 持续 的 精进 ， 避 免 算法 


同时 ，UQSMA 在 保证 精度 的 情况 下 ， 
间 的 方差 更 4 


局 

相 较 于 原始 SMA, 接近 半数 的 测试 函数 结果 分 布 无 明显 差异 ，“” 样 性 、 无 迹 西格玛 点 引入 了 新 的 随机 性 ， 组 合 后 算法 得 益 于 
依 
的 


E 构 原 种 群 ， 保 留 了 原始 SMA 的 搜索 特征 ， 但 在 。” 定性。 直观 来 看 ， 本 文 改进 策略 


反 向 重 构 的 策略 


[hl 


EE 复 计 算 


\， 稳 定性 更 强 ， 说 明 反 向 学 习 扩 大 搜索 空间 的 
特点 有 助 于 保持 算法 的 稳定 ， 即 


反 向 过 程 增加 了 种 群 的 多 


有 优 于 对 比 算 法 的 稳 
组 合 对 防止 算法 停滞 和 加 


强 算 法 稳定 性 有 积极 影响 。 
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录用 定稿 刘 字 准 ， 
为 了 综合 、 简 明 地 评估 UQSMA 与 其 他 算法 对 比 的 竞争 
性 ， 使 用 平均 绝对 误差 MAE 对 上 述 算法 进行 排序 P3。 表 4 
展示 了 所 有 算法 基于 29 个 CEC2017 基准 测试 函数 的 MAE。 
MAE 计算 公式 如 下 : 
和 NF 
MAE- EZ/ | (18) 


只 中 


NF 为 函数 个 数 ，f 为 算法 在 第 i 个 函数 上 求解 结果 的 平 
均值 ， 户 为 第 i 个 函数 的 理论 最 优 值 。 


表 4 MAE 算法 排名 
Tab.4 MAE algorithm ranking 


算法 MAE 


1.80E+05 
2.25E+05 
2.60E+05 
3.12E+05 
1.76E+07 
2.15E+08 
1.10E+09 


UQSMA 
AFSMA 
SMA 
DOLSMA 
HHO 
WOA 
EGSWOA 


值 明 


和 鲁 棒 性 。 


比分 析 不 完全 算法 L291， 
30 将 UQSMA 与 其 不 完全 算法 无 迹 西 格 ] 
法 USMA、 单 纯 使 
单纯 使 
5 展示 了 更 加 综合 和 简 
UQSMA 算法 排名 第 一 
SMA 算法 ,USMA 和 QRSMA 排名 第 三 第 四 日 


算法 


显 小 于 其 他 算法 ， 更 加 综合 


表 4 可 知 ,，UQSMA 排名 第 一 ,并 且 UQSMA 的 MAE 
的 表明 了 该 算法 的 有 效 性 和 


为 了 验证 UQSMA 算法 改进 的 


使 用 CEC2017 基准 测 


实验 二 


H 


有效 忆 
试 函 数 的 维 


E， 需 要 对 
度 


妃 点 弓 


] 拟 反射 学 习 改 进 的 颖 菌 算法 QRSMA 进 
简洁 的 MAE 排名 结果 。 


导 的 秋 菌 算 
j 拟 反问 学 习 改 进 的 黏 菌 算法 QOSMA、 
行 对 比 ， 
表 5 可 见 ， 


表 


，QOSMA 算法 排名 第 二 


优 于 原始 


劣 于 
名 在 单一 使 用 时 


， 说 明了 无 迹 西格玛 点 策 


原始 SMA 


于 引入 新 的 


随机 性 会 导致 算法 精度 下 降 ， 固 所 


破坏 


定 持续 


算法 的 搜索 降低 精度 。 但 所 


名 的 组 合 使 得 算法 表现 最 优 ， 一 方 
的 反 向 重 
的 稳定 性 ， 另 一 方面 无 迹 西 格 


日 | 综 口 


的 组 合 
率 和 


的 拟 反 向 过 程 同样 会 
反射 和 拟 反 向 与 无 迹 西 格 玛 
的 反 向 过 程 和 
FE 西格玛 点 带 来 的 随机 性 ， 


提高 了 


四 点 引入 的 新 搜索 行为 可 


助 算 法 精进 搜索 细致 性 ， 

所 到 精进 广 向 的 盲目 性 ， 防 

， 既 规避 了 单独 策略 使 

使 得 策略 发 挥 出 更 好 的 效果 ， 

稳定 性 的 提升 。 
表 5 MAE 不 完全 算法 排名 

Tab.5 MAE incomplete algorithm ranking 


Hr 


时 的 缺陷 ， 


有 利于 算法 持 


弥补 反 向 学 习 扩 大 搜 
上 算法 停滞 。 多 种 策略 的 


有 索 空 间 后 


同时 策略 
续 搜 索 效 


算法 MAE 


RANK 


UQSMA 1.80E+05 
QOSMA 2.07E+05 
QRSMA 3.45E+05 
USMA 3.64E+05 


一 


果 ， 
要 考 
带 加 


实验 三 : 为 了 验证 UQSMA 算法 在 实际 工程 应 
使 用 UQSMA 求解 轿车 侧面 碰撞 优化 问题 P7， 
虑 了 B 柱 内 径 、B 柱 加 固 、 、 横 梁 

固 、 车 顶 纵 梁 的 厚度 、 护栏 高 立 


i 
pe 
地 
HH 
Hr 
Et 


相 较 


杂 。 


于 其 他 经 工程 优化 问题 该 问题 维 


sa 
度 更 高 、 约 束 更 为 复 


中 的 效 
a 
、 门 


其 目标 函 数 、 约束 条 件 
f(X)=1.98+4.90x, +6.67x, + 6.98x 
十 4.01x4 +1.78xs + 2.73xy 


Subjected to: 


gl(X)=1.16—0.3717x,x4 —0.00931x,mo 
—0.484x,x +0.01343xexio -1<0 


和 决策 变量 取 值 范围 如 下 所 示 。 
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g2(X)=0.261—0.0159xx, —0.188xxs —0.019xx, 
+0.0144x3xs +0.0008757xsxi + 0.08045xcx 
+0.00139xsxi +0.00001575xiox1 -0.32<0 


g3(X)=0.214+0.00817x; —0.131xx, —0.0704xx%, 
+0.03099xxs —0.018xx; + 0.0208x3x8 
+0.121xsx, —0.00364xsxs + 0.0007715xsxio 
一 0.000535xexo+0.00121xx -0.32<0 


g4(X)=0.074—0.061x, — 0.163xsxs + 0.001232xsxio 
—0.166xx%, +0.227x32 -0.32<0 


85(X)=28.98+3.8182 —4.2xx, +0.0207xsxo 
十 6.63x6xo —7.7Xxs +0.32x%Xm0 —32<0 


g6(X)=33.86+2.95x; +0.1792x0 —5.057xx, —11.0x,xs 
—0.0215xsxo —9.98xxs + 22.0xeX0 —32<0 
g7(X)=46.36—9.9x, —12.9xxs +0.1107xx0 —32<0 
g8(X)=4.72—0.5x4 —0.19x,xs —0.0122x4xo 
上 +0.009325xoxo +0.000191x3 -4<0 


g9(X)=10.58—0.674xx, —1.95x,xs + 0.02054xsxo 
一 0.0198x4xio +0.028xexo —9.9<0 


g10(X)=16.45—0.489xsx; — 0.843xsxs + 0.0432xoxio 
—0.0556x,x —0.000786x2? -15.7<0 


where,0.5< A —% <1.5, 
0.192< x —% <0.345, 
—30<xo—X1 <30 (19) 
实验 中 种 群 规模 设置 为 30， 最 大 达 代 次 数 设 置 为 500， 
并 独立 运行 30 次 计算 得 出 统计 结果 。 对 比 用 算法 及 参数 设 
置 与 实验 一 相同 ， 由 于 边界 控制 方式 不 同 EGSWOA 算法 没 
有 使 用 。 求 解 结果 如 表 6 所 示 ， 观 察 可 知 ，UQSMA 算法 虽 
在 重复 计算 中 最 优 值 略 逊 于 原始 SMA 算法 及 其 改进 算法 ， 
但 多 轮 重 复 计 算 的 均值 优 于 所 有 对 比 用 算法 ， 且 方差 较 小 ， 


说 明了 UQSMA 在 实际 的 


化 问题 中 稳定 性 
表 6 轿车 侧面 磁 撞 问题 优化 结果 


和 重 棒 性 更 好 。 


Tab.6 The optimal result for car side impact design problem 
UQSMA AFSMA DOLSMA SMA HHO WOA 
决策 变量 
xl 0.500000 0.500000 0.500540 0.500000 0.500000 0.500000 
x2 1.125390 1.113507 1.121450 1.111922 1.136200 1.157687 
x3 0.500000 0.500000 0.500646 0.500000 0.500000 0.500000 
x4 1.288953 1.307631 1.293401 1.310424 1.499909 1.499801 
x5 0.500000 0.500000 0.500059 0.500000 0.500000 0.500000 
x6 1.499976 1.500000 1.500000 1.500000 1.499909 1.499801 
x7 0.500079 0.500000 0.500000 0.500000 0.500000 0.547264 
x8 0.344995 0.344933 0.345000 0.345000 0.343106 0.192000 
x9 0.330035 0.207846 0.193399 0.305621 0.192000 0.192000 
x10 -17.962397 -20.113941 -18.601569 -20.370221 -30.000000 -30.000000 
xl1 0.840285 -0.000497 -0.117660 -0.550441 -6.681273 5.181629 
优化 结果 
最 优 值 ”22.850268 22.845688 22.848868 22.846324 23.768084 24.040004 
平均 值 ”23.046099 23.211005 23.061354 23.096532 24.901598 27.681259 
标准 差 0.269430 0.355543 0.250977 0.276099 0.792983 1.701868 
约束 值 
gl(X) -0.632695 -0.628194 -0.606523 -0.655046 -0.806886 -0.812774 
g2(X) -0.075343 -0.090981 -0.092269 -0.079360 -0.097477 -0.085725 
g3(X) -0.065285 -0.068270 -0.069034 -0.065982 -0.066378 -0.054259 
g4(X) -0.093732 -0.090221 -0.084601 -0.099203 -0.084642 -0.063955 
g5(X) -3.603556 -4.256299 -4.214778 -3.937901 -5.062355 -4.595731 
g6(X) -6.023978 -7.237602 -7.144709 -6.529494 -9.142410 -7.328384 
g7(X) -0.000799 -0.001838 -0.000935 -0.000773 -0.761912 0.000000 
g8(X) -0.000036 -0.000063 -0.000039 -0.000053 0.000000 -0.005395 
g9(X) -0.936811 -0.974801 -0.948969 -0.978981 -1.140159 -0.820629 
g10(X) -0.276582 -0.185096 -0.158891 -0.264329 -0.217056 -0.341223 


刘 宇 冷 ， 


4 ”结束 语 


针对 SMA 算法 存在 的 不 足 ， 本 文 使 用 了 拟 反 向 和 拟 反 
射 学 习 过 程 组 成 的 综合 反 向 和 无 迹 西 格 玛 点 引导 两 个 策略 对 
原始 SMA 进行 改进 ， 提 出 了 无 迹 西 格 玛 点 引导 的 拟 反 

菌 算 法 UQSMA, 并 将 该 算法 应 用 于 CEC2017 测试 函数 进行 
求解 。 改 进 的 反 向 学 习 框架 将 不 同 的 反 向 生成 过 程 按照 特性 
融合 于 黏 菌 算法 的 勘探 与 开采 过 程 中 ， 扩 大 了 原始 秋 菌 算法 
的 搜索 空间 ， 加 速 了 勘探 阶段 的 收敛 速度 ， 提 高 了 算法 的 稳 
定 
使 


过 


EE 性 ,对 不 同 特点 的 测试 函数 具有 普遍 适用 性 ;迭代 过 程 中 ， 
无 迹 西 格 玛 点 引导 黏 菌 算法 种 群 的 基本 移动 方式 ， 使 其 
移动 方向 更 加 逼近 理论 最 优 解 ， 避 免 算法 的 停 清 。 此 外 ， 为 
了 多 维度 综合 地 评估 改进 算法 的 有 效 性 和 同 其 他 算法 相 比 的 竞 
争 力 ， 本 文 使 用 了 Wilcoxon 秩 和 检验 、MAE 等 指标 ， 同 时 对 
轿车 侧面 碰撞 的 实际 工程 应 用 问题 进行 了 求解 和 比较 ， 多 角度 
地 证 明了 本 文 所 提出 算法 的 求解 精度 、 和 鲁 棒 性 和 实用 性 。 
在 未 来 的 研究 中 ， 和 希望 对 黏 菌 算 法 进行 进一步 的 改进 和 
是 高 该 算法 针对 更 多 样 、 更 复杂 问题 的 适用 性 
其 应 用 于 机 器 学 习 模 型 参数 优化 和 投资 组 合 问 题 求解 等 领域 ， 


步 落实 群 智能 算法 的 实际 应 用 。 
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